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The crystallization process from a solution begins with nucleation, which determines the structure and size 
of the resulting crystals. Further understanding of multi-pathway crystallizations from solution through 
two-step nucleation mechanisms is needed. This study uses density functional theory to probe the 
thermodynamic properties of alumina clusters at high temperature and reveals the thermodynamic 
relationship between these clusters and the saturation levels of dissolved oxygen and aluminum in an Fe-O- 
Al melt. Based on the thermodynamics of cluster formation and the experimental evidence for both excess 
oxygen in the Fe-O-Al melt and for alumina with a polycrystalline structure in solidified iron, we 
demonstrate that the appearance of various types of clusters that depends on the saturation ratio determines 
the nucleation steps that lead to the various crystallization pathways. Such mechanisms may also be 
important in nucleation and crystallization from solution. 

Nucleation and crystallization from solution are fundamental to many industrial processes. Crystallization 
begins with nucleation, which plays a central role in determining the structure, shape and size distri- 
bution of the resulting crystals. Considering the shortcomings of classical nucleation theory (CNT), a 
mechanism for two-step nucleation (TSN) from solution 1 " 5 that includes multi-potential mechanisms 3 ' 6 ' 7 was 
recently proposed and supported by both experimental evidence of stable clusters 1 ' 8 " 15 and by molecular simu- 
lation evidence of the clusters coming together 6,716 " 20 during nucleation. During TSN, sufficiently sized stable (or 
metastable) clusters of solute molecules form first, followed by the reorganization of these clusters into an ordered 
structure 3 during crystallization. The formation of liquid-like or amorphous precursors and even crystal-like 
configurations that differ from the final crystal are essential steps during nucleation 5,21 " 23 . However, the relation- 
ships between the various crystallization pathways and nucleation steps, which are important in controlling the 
structures and sizes of the resulting crystals, are still unclear and need to be studied and supported by experi- 
mental evidence. 

Chemical reactions between dissolved elements in liquid iron that form condensed compounds are very 
common during the melting process. For example, a reaction occurs in which the aluminum (Al) and oxygen 
(O) dissolved in an Fe-O-Al melt form alumina products, usually a-Al 2 0 3 , to decrease the O content and to 
remove as much a- A1 2 0 3 from the melt as possible. The stoichiometric reaction equation is 2A1 + 30 = a-Al 2 0 3 , 
and the change in the standard Gibbs energy is 24,25 A r G m 0 = — 1202.00 + 0.386T, KJ-mol -1 . The process, simply 
described by CNT, is a heterogeneous reaction that results in the nucleation and crystallization of a-Al 2 0 3 . A 
supersaturated Fe-O-Al melt that contains O and Al is required to form a-Al 2 0 3 . The logarithm of the critical 
supersaturation ratio (log S*) was reported to be 3.5 using the electrochemical method 26 , and this result was 
supported by CNT simulations 27 . 

The oxygen content measured following the reaction of Al with O in Fe-O-Al melts frequently exceeds the 
equilibrium level with respect to a-Al 2 0 3 26 " 32 ; however, the oxygen present as non- reactive oxygen continues to 
react with Al to form alumina products that are likely to remain in the solid material during cooling. The excess 
oxygen was hypothesized to originate from both the suspension of fine alumina products 29,30 and from the 
supersaturation of alumina-associated compounds 27,33 . Interestingly, however, the alumina products that were 
observed in solid samples prepared using ultra-rapid cooling at a rate of approximately 10 5 K/s are spherical with 
a diameter of a few to 10 nm or have a coral-like or network-like shape with a trunk thickness of approximately 
several hundred nanometers to a few micrometers, and these alumina products all have polycrystalline structures 
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consisting of a-, y- and 5-alumina 34 . How can the excess oxygen and 
Al react to form products such as metastable y- and 5-alumina 
during the ultra- rapid cooling process but retain their free status in 
the melt? Using CNT and Ostwald's step rule, it was concluded that 
these alumina products may form from the liquid alumina that 
appears during the initial stage of the reaction between Al and O, 
but determining whether liquid alumina is present in an Fe-O-Al 
melt is difficult 35 . 

Using TSN, our work calculates the thermodynamics of the nuc- 
leation of alumina from an Fe-O-Al melt and analyzes the relation- 
ships between the thermodynamics and the experimental evidence. 
Our results provide a clearer picture of multi-step nucleation that 
leads to various crystallization pathways from solution by using the 
thermodynamics of the formation of alumina clusters and the experi- 
mental evidence of both excess oxygen in an Fe-O-Al melt and 
alumina with a polycrystalline structure in solidified iron. 

Results 

Thermodynamic Properties. According to TSN, the formation of 
alumina products in an Fe-O-Al melt occurs as follows: During the 
first step, a population of alumina clusters, (Al 2 0 3 ) n with n = 1-10, 
15 and 30, form when Al reacts with O according to the equation 2A1 
+ 30 = (l/n)(Al 2 0 3 ) n , in which n is the number of A1 2 0 3 units in a 
cluster. The second step is the crystallization of a-Al 2 0 3 from 
(Al 2 0 3 ) n , which occurs according to the equation (1/n) (Al 2 0 3 ) n = 
a-Al 2 0 3 . 

In the current work, the thermodynamic properties of (Al 2 0 3 ) n 
are simulated using density functional theory (DFT). The results 
show that the structures of (Al 2 0 3 ) n (Figure SI a) are in good agree- 
ment with the results of previous studies 36 " 37 . The apparent sizes of 
the (Al 2 0 3 ) n clusters range from 0.410 nm to 1.690 nm (Figure SI b). 
In the temperature range from 0 to 1000 K, the C P (heat capacity), H 
(enthalpy) and S (entropy) of the (Al 2 0 3 ) n clusters increase as the 
temperature increases (Figures S2 a, b, c). The G v (vibrational free 
energy) and G (Gibbs free energy), however, decrease as the temper- 
ature increases (Figures la, b, d). Thus, the change in the thermodyn- 
amic properties of the (Al 2 0 3 ) n clusters with temperature is similar 
to that in the bulk material. Furthermore, the C P , H and S of the 
(Al 2 0 3 ) n clusters increase with increasing n (Figures S2 a, b, c). There 
is a region in which the G v curves intersect in the temperature range 
from 350 to 600 K (Figures la, b). Below the lower limit (350 K) of 
the temperature range, G v increases with increasing n; this trend 
changes above the upper limit (600 K) to G v decreasing with increas- 
ing n. This result indicates that the stability of the (Al 2 0 3 ) n clusters 
changes with variations in the temperature. The zero -point energy, 
E (OK), of the (Al 2 0 3 ) n clusters decreases with increasing 
n (Figure lc). In addition, the thermodynamic properties (C P , 
S, H, Gv) of a-Al 2 0 3 that are simulated using DFT are nearly 
consistent with the experimental data (C P , S) in the temperature 
range from 0 to 1000 K 38 (Figure 2). The E (0 K) of a-Al 2 0 3 is 
-1867993.697 KJ-mol" 1 . 

Thermodynamic Equilibrium. The formation of (Al 2 0 3 ) n clusters is 
characterized using the equilibrium thermodynamics of the Fe-O- 
Al melt. In the second step, that is, (1/n) (Al 2 0 3 ) n = 0t-Al 2 O 3 , A r G m 0 
(S2) is given by A r G m 0 (S2) = G a . A i 2 o3 - (l/n)G (A]2 o3)n- Using G = 
G v + E (0 K), A r G m ® (S2) is calculated as A r G m 0 (S2) = (G v , a -Ai203 
+ E (0 K) a . A12Q3 ) - (1/n) (G v , ( Ai203)n + E (0 K) (A1203)n ). The fitted 
values of A r G m 0 (S2) with respect to temperature in the range of 0- 
1000 K are approximately linear (Table S2, right column); thus, the 
A r G m ° (S2) values in the range of 1850-2000 K are the same as the 
fitted A r G m 0 (S2) values in the range of 0-1000 K (Figure S4 b). In 
the first step, that is, 2A1 + 30 = (l/n)(Al 2 0 3 ) n , A r G m 0 (SI) is given 
by A r G m 0 (SI) = A r G m 0 - A r G m 0 (S2). The A r G m 0 (SI) values in the 
range of 1850-2000 K (which is the steelmaking temperature) are 
also given (Table S2, middle column; Figure S4 a). The results show 



that the A r G m 0 (SI) and A r G m 0 (S2) values are all negative in the 
range of 1850-2000 K; thus, the processes of cluster formation and 
a-Al 2 0 3 crystallization can occur spontaneously. Furthermore, the 
formation of (Al 2 0 3 ) n clusters is more likely than the crystallization 
of 0t-Al 2 O 3 at a constant temperature of 1873 K (Figure S5) because 
the interfacial energy of the melt-liquid interface is considerably 
reduced relative to that of the melt- crystal interface, which results 
in a lower thermodynamic barrier to the formation of (Al 2 0 3 ) n than 
to the crystallization of a-Al 2 0 3 . 

Discussion 

Thermodynamic Behavior of Alumina Clusters. To prove the 
existence of (Al 2 0 3 ) n clusters and to elucidate the nucleation and 
crystallization pathways in an Fe-O-Al melt, the relationship 
between the experimental evidence and the thermodynamics is 
studied. Using the equation 2A1 + 30 = (1/n) (Al 2 0 3 ) n , the 
constant K' can be written as K' = 1/K = (a A i 2 a 0 3 )/a c = [f A i 
(mass% Al) eq ] 2 • [f Q (mass% 0) eq ] 3 /a c , in which K is the equili- 
brium constant; a Ab a Q and a c are the activities of Al, O, and 
(Al 2 0 3 ) n in the melt, respectively; f A i and f Q are the activity 
coefficients of Al and O, respectively; and (mass% Al) eq and 
(mass% 0) eq are the mass fractions (%) of Al and O, respectively. 
In the standard state, the activity of (Al 2 0 3 ) n is 1. The value of f A i 3 • 
f Q 2 is approximately 1 in a dilute Fe-O-Al melt. Thus, K' can be 
written as K' = (mass% Al) 2 eq • (mass% 0) 3 eq . The equilibrium 
relationships between log [mass% Al] and log [mass% O] for the 
formation of (Al 2 0 3 ) n at 1873 K in the present study are plotted as 
a group of colored solid curves (c Q c ) (Figure 3). The measured excess 
oxygen levels 28 ' 30,3 1,39,40 are plotted as a scatter plot using symbols with 
various shapes. The critical value of the supersaturation ratio (S*) for 
the nucleation of a-Al 2 0 3 was measured to be log S* = 3.5 using an 
electrochemical method 26 and was simulated as the log S = 3.5 curve 
and c Q cr value using CNT 27 ; these values are plotted as a dotted curve 
and as a long dashed curve and solid curve, respectively (Figure 3, 
Figure S7). The relationship between Al and O in chemical 
equilibrium with a-Al 2 0 3 and with its associated compound are 
plotted as a dotted curve (c 0 a A1203 > inclined line) and as a short 
dashed curve with a large circle (the equilibrium curve), respec- 
tively 33 (Figure 3). 

The c 0 c curve of (Al 2 0 3 )i is located considerably higher than both 
the c Q cr curve and the experimental data of the excess oxygen from 
various studies. This result indicates that (Al 2 0 3 )i as a single mole- 
cule in an Fe-O-Al melt is unstable and easily transforms into other 
alumina clusters, (Al 2 0 3 ) n with n > 1. The c Q c curves of (Al 2 0 3 ) n 
with n > 1 and most of the measured values of the excess oxygen 
content are located in the region between the c Q cr and Co a ~ A12 ° 3 curves. 
This result indicates that (Al 2 0 3 ) n (n > 1) is present below the 
supersaturation level of the critical nucleus, and the presence of 
excess oxygen means that oxygen is in equilibrium with (Al 2 0 3 ) n 
(n > 1) in the Fe-O-Al melt. 

According to the experimental data of the excess oxygen with the 
c Q c curves of (Al 2 0 3 ) n (n > 1), there are three types of alumina 
clusters in equilibrium with the excess oxygen in the Fe-O-Al melt 
(Figure 3). (I) Alumina clusters. In the cluster region of the c Q c 
curves, the stability of the clusters increases with a decrease in the 
excess oxygen level when the aluminum content (K' decrease) and 
temperature are constant; thus, the clusters in this region are meta- 
stable and change with a variation in the excess oxygen level or 
temperature. (II) Alumina cluster aggregates. In the region between 
the c Q cr and c Q c curves of the most stable cluster, (A1 2 0 3 ) 8 , the clus- 
ters and the critical nucleus coexist. The alumina clusters in the Fe- 
O-Al melt, similar to clusters in a low-temperature solution 3,5,8 , can 
effortlessly coalesce to form aggregates of various shapes and sizes. 
Thus, the excess oxygen in this region should be in equilibrium with 
alumina cluster aggregates that are smaller than the a-Al 2 0 3 critical 
nucleus. Based on this finding, it is concluded that the alumina 



SCIENTIFIC REPORTS | 4:5082 | DOI: 1 0.1 038/srep05082 



2 




Figure 1 | The thermodynamic properties of (Al 2 0 3 ) n with n= 1-10, 15 and 30 at 0-1000 K. (a) The values of the vibrational free energy, G v > of 
(Al 2 0 3 ) n are plotted as a group of colored solid curves, (b) A magnified section of (a) in the region of the intersection of the G v values of (Al 2 0 3 ) n . The 
temperature range of the region of the intersection is approximately 350 K to 600 K. (c) The zero-point energy at a temperature of 0 K, E (0 K). (d) The 
Gibbs free energy, G, where G = G v + E (0 K). 



cluster aggregates are unstable and spontaneously appear depending 
on the degree of saturation of O and Al in the melt. Thus, the sizes of 
the alumina products, which have various shapes, are considerably 
larger than the maximum size of a cluster [( Al 2 O 3 ) 30 , 1 .690 nm] . (Ill) 
Crystal-like clusters. In the region between the c Q c curve of the most 
stable cluster and the c 0 a A1203 curve, the clusters and the 0t-Al 2 O 3 
coexist. Therefore, the excess oxygen in this region should be in 
equilibrium with crystal-like clusters that are more stable than the 
clusters in the melt, similar to the crystal-like clusters in a low- 
temperature solution 3 5 . These results indicate that various species 
of clusters, including alumina clusters, cluster aggregates and crystal- 
like clusters, could be present and in equilibrium with the excess 
oxygen in an Fe-O-Al melt. Thus, we conclude that "the excess 
oxygen" is not really excess oxygen but rather the oxygen in equilib- 
rium with various alumina clusters (hereafter referred to as the equi- 
librium oxygen with clusters) in the melt. 



Multi-step Nucleation Leading to Various Crystallization Path- 
ways. Based on Ostwald's step rule, clusters of various species form a 
series of transitional phases in the process of Al reacting with O. If the 
saturation ratios of Al and O are higher than that of the formation of 
critical nuclei (log S* cr = 3.5, 1873 K, A r G m 0 = -354.169, 
KJ-mol" 1 ), as with Si, there is a multi-step nucleation process in 
the Fe-O-Al melt in which Al and O first react to form (Al 2 0 3 ) n , 
and then, the clusters coalesce to form aggregates; eventually, the 
aggregates crystallize into the critical nucleus, which is followed by 
the immediate and spontaneous growth of a-Al 2 0 3 (1873 K, A r G m 0 
= —479.688 KJ-mol -1 ) during the reaction (Figure 4). However, if 
the saturation ratio is lower than S* cr and higher than that of 
aggregate formation (S* ag ), as with S 2 , the cluster aggregates 
cannot crystallize into the critical nucleus in the melt. In this case, 
the crystallization of 0t-Al 2 O 3 cannot progress further at the 
steelmaking temperature; thus, equilibrium oxygen with clusters is 
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Figure 2 | The thermodynamic properties of the a-Al 2 0 3 crystal. The 

experimental values 38 of C P and S are plotted as blue and green triangles, 
respectively. The values of C P , S, H and G v that were simulated using DFT 
are plotted as colored solid curves. 

present in the post- reaction melt. If the saturation degree is lower 
than S* ag and higher than that of cluster formation (S* d , 1873 K, the 
most stable cluster, (A1 2 0 3 ) 8 , A r G m ® = -381.056 KJ-mol" 1 ), as with 
S 3 , the clusters cannot coalesce and form aggregates in the melt. If the 
degree of saturation is lower than S* d , as with S 4 , the Al and O cannot 
react to form (Al 2 0 3 ) n and there are only crystal-like clusters in the 
Fe-O-Al melt (Figure 4). We conclude that during cooling, various 
species of clusters in the post-reaction Fe-O-Al melt can only 
crystallize directly into alumina with different shapes, sizes and 



polycrystalline structures that consist of a-, y- and 5-Al 2 0 3 
because the equilibrium thermodynamics change with the decrease 
in temperature (Figure S10). In other words, the morphologies and 
structures of the alumina products in solid steel could strongly 
depend on the various species of clusters in the melt. 

Because of these findings, we argue that stable clusters could be 
present in an Fe-O-Al melt. In a melt with a high saturation ratio, 
alumina clusters can coalesce and form aggregates that then change 
into critical nuclei, leading to the crystallization of a-Al 2 0 3 during 
the reaction of aluminum with oxygen. In a post-reaction Fe-O-Al 
melt, various species of alumina clusters are present depending on 
the saturation ratio, and they cannot crystallize into the critical nuclei 
of a-Al 2 0 3 but rather crystallize directly into various structural alu- 
mina products (y-, 5-Al 2 0 3 ) as the temperature decreases during 
the cooling stage. A greater understanding of the multi-step nuc- 
leation that leads to the various crystallization pathways may lead 
to the design of auxiliaries that can control the rate of nucleation and 
avoid the formation of undesired solid forms. 

Methods 

All the simulations for calculating the thermodynamic properties in the present work 
are performed using DFT with the molecular orbital theory computational program 
Dmol 3 . The initial structures are established using the Visualizer module of MS 6.0. 
The geometry optimization is performed using the BFGS (the abbreviation corre- 
sponds to the first letters of the names of the following researchers: Broyden, Fletcher, 
Goldfarb, and Shanno) method with a quasi-Newton algorithm. The hybrid density 
functional BLYP (the abbreviation corresponds to the first letters of the names of the 
following researchers: Becke, Lee, Yang and Parr) method using a generalized gra- 
dient approximation (GGA) is used as the exchange-correlation potential function. 
The thermodynamic properties of various structures are calculated using the atomic 
harmonic vibrational frequency based on statistical thermodynamics. The initial 
structures of the alumina clusters, (Al 2 0 3 ) n with n = 1-10, 15 and 30, are established 
using Visualizer. The precision of BFGS is set as follows: energy < 2.0 X 10~ 5 Ha, 
tension < 0.004 Ha/A, and shift < 0.005 A. The self- consistent field (SCF) method is 
used with the precision of the total energy and charge density set at 1 X 10" 5 Ha, and 
the thermal smearing effect is used at a precision of 0.005 Ha. The cut-off radius of the 
DNP basis set of the d orbital is 3.5 A. Electrons outside the atomic nucleus are 
handled using the effective core potentials (ECP) method. The initial structure of the 
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Figure 3 | The relationship between log (mass% Al) and log (mass% O) at 1873 K. The equilibrium relationships between (Al 2 0 3 ) n and Al and O in an 
Fe-O-Al melt are plotted as a group of colored solid curves (c 0 c ). The measured values of excess OX yg en 28 > 30 ' 31 > 39 > 40 are plotted as a scatter plot using symbols 
with various shapes. The critical value of the supersaturation ratio (S*) of oc-Al 2 0 3 nucleation was measured to be log S* = 3.5 using an electrochemical 
method 26 and was simulated as the c Q cr value using CNT 27 ; these values are plotted as a dotted curve and as a long dashed curve and solid curve, 
respectively. The relationships of Al and O in chemical equilibrium with a-Al 2 0 3 and with its associated compound are plotted as a dotted curve (c 0 aA1203 , 
inclined line) and as a short dashed curve with a large circle (the equilibrium curve), respectively 33 . 
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Figure 4 | A schematic illustration of the nucleation and crystallization from the melt. S* cr is the critical degree of saturation at which the critical nuclei 
form. S* ag is the critical degree of saturation at which the cluster aggregates form. S* c i is the critical degree of saturation at which the alumina clusters 
form. The melt in the SI state can form critical nuclei and crystallize into alumina products, that is, a-Al 2 0 3 in an Fe-O-Al melt, and the multi-step 
nucleation process is shown as a solid red curve. However, S 2 , S 3 and S 4 can only form alumina products that are dependent on the crystallization process 
during cooling of the post-reaction melt, and the multi-pathways to crystallization are shown as solid black curves. 



a-Al 2 0 3 crystal is adopted from the MS 6.0 structural database. The precision of the 
BFGS algorithm is set as follows: energy < 2.0 X 10" 5 Ha, tension < 0.004 Ha/ A, and 
shift < 0.005 A. The Brillouin zone integral is calculated using the Monkhorst_Pack 
method with a grid size of 3 X 3 X 2 k-points. The precision of the total energy and 
charge density is set as 1 X 10" 5 Ha, and the thermal smearing effect is used at a 
precision of 0.055 Ha. The cut-off radius of the DNP basis set of the d orbital is 3.5 A. 
Electrons outside the atomic nucleus are handled using the ECP method. 
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